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Nomenclature 

A cross-section area, cm- 

C relative atom concentration 

En nitrogen dissociation, MJ/cm^ 

Hj total enthalpy, MJ/kg 

hp> enthalpy of formation, MJ/kg 

I current, amps 

k constant 

k w reaction rate constant, cm/s 

L e Lewis number 

M molecular weight, gm/mole 

m mass loss rate, gm/s 

mf mass flow rate, gm/s 

P pressure, atm 

Pr Prandtl number 

q heat flux, W/cm 2 

R radius, cm 

91 gas constant, atm cnvVgm mole K 

9? e 2 Reynolds number 

Sc Schmidt number 

s arc length, cm 

T temperature, K 

U velocity, cm/s 

V voltage, volts 

a mass fraction 

y a absolute recombination coefficient 


Yi apparent recombination coefficient 

e emittance 

p viscosity, N-s/m 2 

p density, gm/cm^ 

o Stefan-Boltzmann constant 

Subscripts 


A 

air 

abs 

absolute 

ch 

chamber 

e 

boundary layer edge 

eff 

effective 

F 

flat-faced cylinder 

f 

frozen 

H 

hemisphere 

i 

chemical species 

N 

nitrogen 

0 

oxygen 

0 

stagnation point 

th 

total hemispherical 

w 

wall 

oo 

free stream 

2 

behind bow shock wave 

Superscript 

* 

sonic point in nozzle throat 
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Summary 

Surface properties have been obtained on several classes 
of thermal protection systems (TPS) using data from both 
side-arm-reactor and arc -jet facilities. Thermochemical 
stability, optical properties, and coefficients for atom 
recombination were determined for candidate TPS 
proposed for single-stage-to-orbit vehicles. The systems 
included rigid fibrous insulations, blankets, reinforced 
carbon carbon, and metals. Test techniques, theories used 
to define arc -jet and side-arm-reactor flow, and material 
surface properties are described. Total hemispherical 
emittance and atom recombination coefficients for each 
candidate TPS are summarized in the form of polynomial 
and Arrhenius expressions. 

Introduction 

This report describes the research conducted as part of the 
Access-to-Space phase I study to develop a database on 
surface properties of candidate thermal protection systems 
(TPS) for a single-stage-to-orbit (SSTO) vehicle (ref. 1). 
Candidate systems include rigid fibrous insulations, 
flexible blankets, reinforced carbon carbon, and metallics. 

During atmospheric entry, at hypersonic speed, the high 
temperature air between the bow shock wave and the TPS 
surface will be partially or fully dissociated into atoms. 
Therefore, the heat transferred to the surface of the TPS 
will consist of chemical as well as sensible energy. The 
chemical energy transferred (catalytic efficiency) to the 
surface of the TPS can strongly influence the heat transfer 
rate to the vehicle. For example, the Orbiter surface 
catalysis experiment demonstrated a 40% eflect on the 
heat transfer rate during Earth entry (refs. 2-4). There- 
fore, to accurately size a TPS for any proposed SSTO 
vehicle, using one or a combination of materials, the 
designer must include the surface catalytic efficiency 
along with other surface properties such as thermo- 
chemical stability and emittance (refs. 5-9). 

In this study, thermochemical stability, emittance, and 
atom recombination coefficients for 16 candidate SSTO 
TPS materials are reported. The Laser Accurate Surface- 
catalytic Efficiency Research (LASER) complex, test 


techniques, and theories developed at Ames Research 
Center to obtain these surface properties and calculate 
atom recombination coefficients are reviewed 
(refs. 10-12). 

The author gratefully acknowledges the support of 
Y. K. Chen in obtaining the BLIMPK code solutions, and 
of J. Marschall and J. Pallix in obtaining the coefficients 
from the side-arm reactor. These tasks were supported 
under Eloret Institute contract NCC2-462. Laser-induced 
fluorescence diagnostic measurements were conducted by 
Douglas Bamford under a NASA Small Business 
Research Program contract NAS2- 13469. 

Complex 

Ames Research Center has the unique capability of being 
able to obtain the surface properties, including catalytic 
efficiency, of candidate materials from a single laser 
complex located at Ames (fig. 1 ). The complex consists 
of a dye laser, a side-arm-reactor facility, and an arc-jet 
facility that are interactive. The central location of the 
laser provides for the use of laser-induced fluorescence 
(LIF) diagnostic techniques to determine the relative 
specie concentration levels in both the side-arm reactor 
(fig. 1(a)) and free-arc -jet streams (fig. 1(c)). The 
complex greatly reduces the cost and time required for 
obtaining surface properties of candidate TPS materials 
from room temperature to their upper use temperature. 

Laser Facility 

The dye laser used to excite the atoms in the tube 
containing the sample is located in a room next to the 
side-arm-reactor facility (fig. 1(b)). A schematic of the 
laser system integrated with the side-arm-reactor compo- 
nents is shown in figure 2(a). The system was developed 
by Bamford and Romanovsky (ref. 13). A frequency- 
tripled Nd:YAG laser (Continuum NY81) pumps a 
tunable dye laser (Continuum ND60, Dual Grating 
Option) to excite oxygen or nitrogen atoms in either of 
the facilities. To detect oxygen atoms, Coumarin 460 
laser dye is used to produce tunable radiation in the 452 
nm region. For nitrogen, Stilbene 420 dye is used to 



produce tunable radiation near 422 nm. The beam from 
the dye laser is passed through a doubling crystal of 
beta-barium borate (BBO) to produce the appropriate 
ultraviolet wavelength. The ultraviolet radiation is 
separated from residual blue radiation with a Pellin-Broca 
prism and is reflected off several dielectric mirrors and 
quartz prisms toward the arc-jet or side-arm-reactor test 
cell. The beam is used to excite the oxygen or nitrogen 
atoms in the test region via a two-photon transition at 
either 226 or 21 1 nm, respectively. Fluorescence by the 
atoms is detected using one or more gated photomultiplier 
tubes in conjunction with the appropriate narrow band 
filter for wavelengths at 845 nm (oxygen) and 869 nm 
(nitrogen). Signals are collected using gated integrators 
and are stored on a computer. 

Side-Arm-Reactor Facility 

The side-arm-reactor facility (Fig. 1(c)) essentially 
consists of a gas delivery and pressure control system 
coupled to a tubular reactor made up of a main arm and a 
dead-end side arm. High purity gases are fed into the 
reactor through a metering valve and the gas flow is 
measured using Tylan mass How meters. The system is 
evacuated using turbo-molecular and roughing pumps. 

The system pressure is monitored using Baratron-type 
capacitance manometer gauges, and the target pressure is 
obtained by adjusting the gas inflow. 

A microwave-discharge cavity located just upstream of 
the side-arm section was used to dissociate the test gas 
into atoms. A hinged heavy-duty electric furnace (68 cm 
long), positioned around the reactor section (side arm) of 
the facility, was used to raise the temperature of the test 
samples. The heater was located several centimeters away 
from the Tee junction formed by the main line and the 
side-arm section. The wall temperature within the side 
arm can be varied from 300 to 1 300 K. A typical test 
specimen (1 to 27 cm long) is positioned at the center of 
the heater. Either a thermocouple probe (ref. 12) or LIF 
diagnostics can be used to obtain the data necessary for 
calculating the atom recombination coefficients (ref. 1 1). 

Earlier coefficients were calculated using a thermocouple 
probe (consisting of a platinum platinum/ 13% rhodium 
thermocouple and coated with either manganese or silver) 
to measure the temperature distribution along the center- 
line of the tube containing the test specimen. The repro- 
ducibility of the thermocouple measurements was within 
1%. The present method uses LIF measurements to obtain 
the distribution of the specie concentration through this 
area of the tube. The atoms are excited by the laser energy 
focused through the end of the side arm, and the radiation 
from the excited atoms is measured using photomultiplier 
tubes which are aligned with the ports in the side wall of 


the heater. The reduction in atom concentration along 
the tube is used to determine the atom recombination 
coefficient. 

Arc-Jet Facility 

The typical measuring equipment used to obtain data 
from the aerothermodynamic heating facility (AHF) at 
Ames Research Center is shown in figure 2(b). The AHF 
uses a constricted arc heater to provide high-enthalpy 
dissociated hypersonic flow over a test model positioned 
downstream of a 16 deg conical nozzle. Either nitrogen 
or air can be easily used as the test gas without altering 
the heater hardware. This permits quick, consecutive 
measurements of heat flux and temperature to be taken 
from a test model during its exposure to either test gas. 
Surface conditions on the test model are varied by 
changing the exit diameter of the nozzle, the reservoir 
pressure, or the electrical power dissipated in the arc 
heater. The geometric area ratio (nozzle exit to throat) of 
the facility can be varied from 64 to 400. Heater pressure 
can be varied from 0.68 atm to roughly 5.5 atm and the 
maximum power dissipation in the heater can be 
increased up to 20 MW. Stagnation point enthalpy was 
determined using a nozzle code in conjunction with 
velocity measurements taken using a LIF diagnostic 
technique (ref. 10). 

Test Articles 

Atom recombination coefficients were obtained by 
conducting tests on samples of the actual candidate TPS 
in both the arc -jet and side-arm-reactor facilities. 

Side-Arm Reactor 

Split tube or short tubular sections of the proposed TPS 
were used in the side-arm reactor to obtain data needed to 
calculate the recombination coefficients (fig 1(c)). The 
thermal control coating was applied to the inside surface 
of split tubes made from materials such as rigid fibrous 
insulation or carbon. Fabric and metallic materials were 
tested as liners inside a quartz tube. Depending on the 
material's relative catalytic efficiency, tube lengths 
ranged from 2.5 cm to as much as 22.86 cm. The inside 
diameter of the tube was roughly 1 .9 cm. 

Arc Jet 

The arc -jet test configurations were designed to provide 
an adiabatic back wall and ensure uniform temperature 
and pressure distributions across the test samples. During 
the arc-jet tests, two different test configurations were 
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used to obtain data (fig. 3): a 5 deg flat-faced cone and a 
15.2 cm diameter flat-faced cylinder. The flat-faced cone 
was made using rigid fibrous insulation (fig. 3(a)). It had 
an 8.3 cm base diameter, a 1 .3 cm corner radius, and a 
thickness of 6.35 cm. A high emittance surface was 
developed by either coating or impregnation of the cone. 

A threaded aluminum mounting ring was bonded into 
the base of each cone so that it could be attached to a 
water-cooled support. Surface thermocouples (platinum/ 
platinum 13% rhodium) were installed near the stagnation 
point of the cone. 

The second configuration used samples cut in the shape 
of a 7.1 1 cm diameter disk. They were positioned at the 
center of the model by being mounted inside a retaining 
ring and backed by a disk of rigid fibrous insulation. The 
retaining ring, 7.62 cm in diameter and 6.25 cm thick, 
was also made from a coated rigid fibrous insulation. 
Finally, platinum/platinum/ 13% rhodium thermocouples 
were installed behind the sample, either in a plug (coated 
with a borosilicate glass) or in the surface of the backing 
material. This model holder resulted in the sample being 
recessed 0.3 cm below the front surface of the cylinder. 
Earlier arc-jet tests showed that the recessed mounting of 
the sample did not affect the surface temperature or heat 
flux relative to a flush-mounted sample (ref. 14). 

Candidate TPS 

Candidate TPS materials used in this study are listed in 
table 1. 

All candidate systems tested using the cone configuration 
(fig. 3(a)) were developed by Ames Research Center. 
These include: 

1 . Reaction-cured glass (RCG) (ref. 15) applied over a 
cone made using a mixture of silica and alumino- 
borosilicate fibers. RCG is basically a fully dense 
borosilicate glass roughly 0.030 cm thick. 

2. Toughened unipiece fibrous insulation (TUFI), 
which consists mainly of a mixture of molybdenum 
disilicide and borosilicate glass (ref. 16). This system 
was impregnated into the surface of the insulation to 
a depth of 0.254 cm before being sintered in the 
furnace. The density of TUFI is roughly half of RCG. 
For this study, TUFI was impregnated into AETB-12, 
a mixture of aluminoborosilicate, silica, and alumna 
fibers. 

3. Silicone-impregnated reusable ceramic ablator 
(SIRCA) made with a high purity silicon impreg- 
nated into a billet of FRCI-12 (a mixture of silica and 
aluminoborosilicate fibers) (ref. 17). Since SIRCA is 
an ablator (active system), its surface properties are 


dependent on the aeroconvective heating 
environment used to form the char. 

Flexible blanket systems (ref. 18) included in the study 
are: 

1 . Tailored advanced blanket insulation (TABI), 
developed by Ames Research Center, made using 
silicon carbide or aluminoborosilicate fibers 
(Nextel 440). 

2. The Nextel 440 system coated with the protective 
ceramic coating (PCC), also developed by Ames 
Research Center. PCC basically consists of Ludox 
and uses silicon tetraboride (SiB6) as a flux and 
emittance agent. 

3. The Nextel 440 system coated with the shuttle grey 
C-9 coating, developed by Rockwell International. 
The C-9 coating was made using tetraethyl 
orthosilicate (TEOS) glass along with a silicon 
carbide (SiC) emittance agent. 

Coated advanced reinforced carbon carbon systems 
were provided by Langley Research Center, ROHR 
Incorporation, and McDonnell Douglas Aerospace 
(MDA). These coatings include: 

1 . C-CAT ACC-4, made by LTV and Carbon Advanced 
Technologies, which basically consists of a carbon 
composite with an SiC conversion coating treated 
with a TEOS glass and covered with a Type 1 sealant 
(sodium silicate-based glass). 

2. LVP, a chemical vapor deposition (CVD) SiC/SiB 
glass containing zirconium oxide, which was made 
using a Loral and Vought Process. 

3. ORCC (RT-42), obtained from ROHR, is a three- 
part system which consists of an initial sealer, an 
oxidation barrier (near stoichiometric mixture of 
B4C + SiC that was applied using CVD), and a 
borosilicate glass overglaze sprayed onto the 
samples. 

4. Silicon carbide-silicon boride coated system, made 
by the Russians using CVD, was obtained from 
MDA. Finally, a sample of an uncoated SiC/SiC 
composite was provided by MDA for the study. 

Surface sheets of the following candidate metallic TPS 
were characterized during this investigation: 

1 . Oxidized and unoxidized samples of Inconel 617 
(nickel-based material). The surface layer on the 
oxidized samples was the result of exposure for 2 hr 
in air at a temperature of 1256 K ( 1 800°F) and a 
pressure of 1 atm. These samples were provided 
by MDA. 
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2. Oxidized samples of MA-956 (an iron-based 
material). 

3. PM- 1000, with a composition similar to Inconel 617, 
was provided by ROHR. 

Experiment 

Porosity Measurements 

Open porosity of the surface on the rigid and flexible 
fibrous insulations is directly proportional to the 
material’s effective surface area; therefore, it directly 
affects the surface catalytic efficiency. The open surface 
porosity of each material was determined using an 
apparatus illustrated in sketch 1. 



Sketch 1 


Using a very low mass-flow rate through the apparatus, 
the expression for conservation of mass, and Bernoulli’s 
equation, the open-surface porosity of a TPS system can 
be calculated from the following equation; 

A etf =Cxm,/D 2 (APxP t ) l/2 (1) 

The calculated open-surface porosity for the silicon 
carbide cloth system was greater than 50%. The porosity 
of the other cloth systems and the TUFI was below 10%. 
Porosity results in increased surface area, thereby 
increasing the surface catalytic efficiency above a value 
for a solid surface. The coefficients calculated from the 
data collected from both the arc jet and side-arm reactor 
are unique because they were obtained using the actual 
proposed TPS materials. Therefore, these apparent 
recombination coefficients include the effect of surface 
area (porosity) for each TPS. Earlier experiments carried 
out on a quartz surface in a side-arm reactor showed that 
the coefficient y, increased linearly with surface area 
(ref. 19): 

Yi = Tabs (2) 

Apparent coefficients establish the relative catalytic 
efficiency between materials and are required in order to 


determine the total heat flux to the vehicle during high- 
speed Earth entry. The correction for surface area (S a ) in 
equation (2) is equal to the ratio of the chemical active to 
the geometric surface areas. It should be made clear that 
many of the recombination coefficients reported in the 
literature and those presented here are not absolute values, 
but apparent values (y,). 

Side-Arm-Reactor Tests 

Tests in the side-arm reactor were conducted at constant 
pressure but at temperatures ranging between room 
temperature and 1250 K. The size of the power supply 
and microwave cavity was chosen in order to provide 
test conditions over a range of pressures from 0. 1 to 
1.0 torr. These experiments were conducted at 0.27 torr 
(3.5 x lO^ 4 atm) using both high-purity molecular oxygen 
and nitrogen test gases. Data taken from this facility to 
determine the coefficients for the candidate systems were 
obtained over surface temperatures ranging from room 
temperature to 1250 K. Data collected from cloth systems 
in earlier studies were limited to much lower temperatures 
than the present study because of outgassing of the test 
specimen and the interaction of the volatiles with the 
thermocouple at this very low pressure (ref. 5). 

To reduce the effect of outgassing on the data obtained 
from the cloth system, it was preconditioned in a vacuum 
furnace at 900 K and 2 mmHg for 24 hr to burn out the 
sizing material before exposure to the low-pressure 
side-arm-reactor environment. Using LIF diagnostic 
techniques and the preconditioned sample increased the 
temperature at which data could be taken from roughly 
500 K to over 1200 K for the cloth systems. 

Arc-Jet Tests 

Tests were also conducted using both nitrogen and air 
streams in the arc -jet facility. The models were exposed 
to the stream for 180 sec to ensure steady-state surface 
conditions during the measurements. Typical test 
conditions are shown in table 2. Arc-jet chamber pressure 
was varied from 0.68 to 4.1 atm and the power dissipation 
in the gas ranged from roughly 1 .0 to 7.2 MW. Finally, a 
60.96 cm diameter nozzle exit was used in this study 
(geometric area ratio between the nozzle exit and throat 
was 256). In addition, the mass flow through the heater 
was calculated using a fixed orifice size between the gas 
supply and the heater manifold. Total enthalpy (bulk) of 
the reservoir was then calculated using the following 
power and mass flow relationship: 

H T =qkIV/(m f ) (3) 
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where 

q = efficiency factor (W. L. Love, private 
communication, 1987) 

k = conversion constant 

Note that the bulk enthalpy was always lower than the 
effective enthalpy determined from the free-stream 
velocity measurements (table 2). 

During each test the radiated energy from the stagnation 
region (roughly a 0.6 cm diameter area on the model) was 
measured using a radiometer with a bandwidth from 0.2 
to 10.2 (I. Stagnation point pressure on the test models, 
located on the centerline and 35 cm downstream from the 
nozzle exit, was inferred from measurements taken with 
a water-cooled hemisphere. The surface pressure was 
measured using a bridge-type transducer with an accuracy 
of ±5% full scale. In addition, pyrometer and thermo- 
couple measurements were made to determine the surface 
temperature of each test sample. Samples were tested 
from one to five times over a range of surface tempera- 
tures from 1200 to 1730 K. Surface pressures ranged from 
0.005 to 0.035 atm and effective enthalpies from 12.8 to 
25 MJ/kg. 

Pre- and posttest data were obtained from each sample 
in the form of photographs, reflectance measurements, 
and X-ray fluorescence analysis. These data included 
(1) room temperature spectral reflectance measurements 
using a BIO-RAD model FTS 40 (wavelength 0.25 to 
2.5 microns) spectrophotometer and a Perkin Elmer 
model 310 (wavelength 2.5 to 20 fi) spectrophotometer 
and (2) elemental chemical analysis of the surface of each 
sample using X-ray fluorescence measurements. 

Finally, LIF diagnostic techniques (ref. 5) were used to 
aid in defining the free-stream properties in the 
hypersonic arc-jet streams. 

Analysis 

Low-Temperature Coefficients (Side-Arm Reactor) 

Low-temperature coefficients were obtained over a 
temperature range from room temperature to 1250 K 
using the side-arm reactor. In the past, coefficients were 
determined using data obtained from a thermocouple 
probe mounted along the centerline of the tube 
(refs. 17-19). The probe temperature increases because 
of atom recombination on its surface and is directly 
proportional to the local atom concentration. The 
temperature of the probe decreases as the probe is moved 
along the centerline of the tube, within the sample. For a 
first-order surface reaction, the temperature difference 


between the probe and sample will decrease exponentially 
down the length of the tube (ref. 17). 

A typical temperature distribution for the low catalytic 
efficiency material is illustrated in sketch 2. 



AT = K xexp(-mL/ R) (4) 

where 

t / 'j 

m =(VmRk^ /2.D12) 

Equation (4) shows a relationship between the 
temperature difference (AT) and the distance along the 
tubular-shaped sample, derived by Smith (ref. 20) and 
Shuler and Laidler (ref. 21 ), to calculate the atom 
recombination coefficient. The constant (K) in equa- 
tion (4) depends on the heat transfer rate, gas transport 
properties, and probe surface catalytic efficiency. The 
recombination coefficient is directly related to the square 
of the slope of the ln(AT) versus L/R curve. Diffusion 
coefficients for atomic oxygen and nitrogen were 
empirically derived from the Chapman-Enskog equation 
using values of collision cross sections calculated by Yun 
and Mason (ref. 22). The thermal properties for these 
calculations were obtained from a thermal model derived 
by Cagliostro (private communication, 1994). 

Currently, the atom recombination coefficients are 
calculated using measurements from LIF diagnostic 
techniques to obtain the specie concentration profile 
along the tube centerline. This technique was first 
demonstrated by Pallix and Copeland (ref. 23) tor a 
relatively low catalytic efficiency material, a long tubular 
sample (L/D > 10). Typical PMT signals and relative 
oxygen atom concentration profiles for both low and 
moderate catalytic efficiency materials are shown in 
figure 4. The signal profile for the low catalytic efficiency 
material (fig. 4(a)) is similar to the temperature profile 
shown in sketch 2. The coefficient can be calculated from 
the diffusion expression in equation (4) by substituting 
the slope of the PMT signal for AT. 
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PMT signals representative of specie concentrations along 
the centerline of the tube for a moderately high catalytic 
efficiency material (silicon carbide/glass-coated cloth) are 
plotted relative to time and distance (fig. 4(b)). For this 
material the specie profile is no longer exponential along 
the centerline of the tube. Note that because of the high 
rale of atom recombination on its surface a much shorter 
sample was required. These data were obtained using a 
three-section diffusion tube (short sample of test material 
(L/D < 1) positioned between two sections of quartz tube) 
which was placed at the center of the heater. Therefore, to 
calculate atom recombination coefficients for moderate 
and high catalytic efficiency materials, the basic diffusion 
equation must be solved using appropriate boundary 
conditions at the tube interfaces (ref. 12). 


d 2 c lac a 2 c 

dr 2 + r3r + 9*2 


= 0 


(5) 


where r and x are dimensionless cylindrical coordinates 
measured in units of tube radius and C = C(r,x)/C 0 is the 
atom concentration normalized by the value at x = 0. This 
approach uses the same assumptions as the earlier method 
using the thermocouple: (1 ) gas phase recombination of 
atoms is negligible, (2) recombination of atoms on the 
wall of the sample is a first-order reaction, and (3) the 
side arm is at uniform temperature. 

The LIF methods are more direct and less intrusive than 
using the thermocouple probe to obtain data for calcu- 
lating the coefficients. A detailed description of the 
solutions for moderate and high efficiency materials is 
given in reference 12. 


High-Temperature Coefficients (Arc Jet) 

Surface heat flux to the stagnation point of a test model 
exposed to arc -jet flow can be calculated using either a 
full Navier-Stokes code such as the General Aerodynamic 
Simulation Program (GASP) (ref. 24) or, if the flow is 
chemically frozen, a theory such as Goulard’s (ref. 25). 
Computational differences in the two solutions are shown 
schematically in figure 5. 

GASP was applied as an axisymmetric nonequilibrium 
air chemistry solution of the Navier-Stokes equations 
using a one-temperature model (fig. 5(a)). The two-factor 
approximate factorization algorithm and Van Leer flux 
vector splitting method (with upwind-biased third-order 
differencing) were incorporated into the solution using a 
global iteration scheme. A five-species (N, O, N 2 , O 2 , 
and NO) chemical reaction model was used in the code 
for air. The reaction rates applied in the solutions are from 
the work of Park (ref. 26). 


In order to calculate the heat to the stagnation point on a 
test model, two GASP solutions are required: one for the 
nozzle flow and one for the flow through the bow shock 
wave and over the model. Therefore, axisymmetric grids 
are required to represent the flow through the nozzle and 
lest section and the flow around the model. The code 
requires inputs of the heater conditions, atom recombina- 
tion coefficients, and surface emittance for the materials 
during the arc -jet test to calculate the heat flux to the 
model. Also, the location of the model from the nozzle 
exit must be known. The formulation of this code is 
presented in detail in the GASP User’s Manual (ref. 24). 

Unlike GASP, Goulard’s theory requires that the flow 
properties in the solution be obtained independently 
using a nozzle code. The general equation for Goulard’s 
theory is: 

q wo =0.66Pr- 2/3 (p2H2) l/2 [(du e /ds) F ] [H e0 -H w ] 

l + (Le 2/3 <|> 0 -l)a 0 h D (Le 2/3 <)> N - l)a N h D 

X (H 0 -H w ) + (H 0 -H w ) 

( 6 ) 


where 


|1 + 


0.47S" 2/3 [2(du e /ds) F *p 2 |i2 


Pw ^wi 


Inputs required in Goulard’s theory (based on a 
hemispherical-shaped configuration) are the gas 
properties from behind the bow shock wave and at the 
stagnation point of the model. These properties include 
the velocity gradient, wall reaction rate constant, and 
enthalpy. When the flow properties and velocity gradient 
are known, the reaction rate constants for a given material 
can be calculated from the equations by using a measured 
stagnation point heat flux and surface temperature as 
inputs. 


SCFC Code 

In this study, a code was developed to calculate atom 
recombination coefficients and the total hemispherical 
emittance for advanced TPS using arc -jet data taken from 
test models in the shape of either a blunt cone (fig. 3(a)) 
or flat-faced cylinder (fig. 3(b)). Surface coefficients are 
calculated assuming frozen chemistry (SCFC) and incor- 
porating Goulard’s theory as part of a nozzle program 
written by Yoshishawa and Katzen (ref. 27). The program 
is based on the assumption that most properties of the gas 
in the test section (downstream of the nozzle), having 
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undergone relaxation of several internal degrees of 
freedom, are approximately the same as if the gas had 
made an instantaneous transition from full equilibrium 
flow to flow in which all internal energy exchange is 
frozen (fig. 5(b)). In one-dimensional flow, this process 
can be characterized by a parameter called the frozen 
Mach number (Mf), which is used to define the state of 
the gas in the free stream. 

In order to determine Mf, the effective area ratio of the 
test volume to nozzle orifice must first be estimated. 
Because impact pressure is relatively insensitive to Mf, 
it can be used along with the calculated bulk enthalpy 
(heater) to estimate the effective area ratio (A/A*). On the 
other hand, free-stream velocity, temperature, and atomic 
specie concentrations (oxygen and nitrogen mass fraction) 
in the flow are very sensitive to Mf. Figure 6 shows that 
the velocity increases by over 800 m/s and the nitrogen 
mass fraction (in both dissociated air and nitrogen 
streams) decreases by as much as 20% (0.2) between the 
total frozen chemistry condition (Mf = 1) and equilibrium 
flow. In addition, for the same reservoir conditions, the 
amount of atomic nitrogen (mass fraction) present in a 
nitrogen stream is higher than for air. 

Free stream- During this study, the flow properties 
(frozen Mach number) were obtained by iterating between 
the enthalpy and velocity, found from the LIF diagnostic 
techniques, until a minimum difference or agreement 
between measured and predicted velocity was reached 
(fig. 7). The resulting total enthalpy (effective enthalpy) 
was consistently 5% to 15% above the calculated bulk 
enthalpy values obtained using equation (3). However, 
the variations in enthalpy over the test conditions used for 
these tests were well within the performance envelope of 
the heater. The properties calculated using the LIF 
diagnostic technique will be referred to in this paper as 
“measured” values. 

LIF-measured free-stream velocities were determined 
using the Doppler shifting of the fluorescence peaks of 
the nitrogen atoms and oxygen atoms in the nitrogen and 
air flows, respectively (fig. 7). Good agreement between 
measured and predicted values from the SCFC code was 
achieved for the velocity in most test cases for nitrogen. 
However, at the two lower test cases in air (P c h = 0.68 
and 1.36 atm), the measured values could not be corre- 
lated any closer than about 5% to the predicted values 
(fig. 7). Further study is required in order to resolve the 
differences between the measured and predicted values 
because of uncertainty in the heater performance at low 
pressures. In addition, the frozen Mach numbers deter- 
mined for the test cases were also substantiated by 
comparing measured and predicted values of nitrogen 
specie concentrations and translational temperatures in 


air (figs. 8 and 9). Atomic nitrogen mass fraction in the 
hypersonic streams was plotted relative to the available 
amount of chemical energy per unit volume of gas 
available to dissociate molecular nitrogen during each 
test. The chemical energy per unit volume was calculated 
using the following expression: 


e n 


|h T -|k/2(U 00 ) 2 



(7) 


Atomic nitrogen concentrations in the free stream were 
plotted as normalized values relative to the amount of 
atomic nitrogen found to exist during test case 1 in air. 
Predicted values from the SCFC code account for the 
differences in the mass flow rates between test cases. In 
general, the predicted and measured values agreed well 
with each other for the air tests. 

Also, measured values of translational temperature were 
calculated using the Gaussian line shapes (produced by a 
two-photon transition) that were recorded during the LIF 
experiments (fig. 9). Included in the figure are values 
assuming equilibrium flow and curves for Mf = 1 .4 and 
Mf = 2.5 using the SCFC code. The SCFC calculation 
assumed that the translational, rotational, and vibrational 
temperatures were equal. The figures show clearly that 
the chemical state of the gas in the free stream during 
these tests correlates well with the predicted values near 
Mf = 2. 

Within the accuracy of the LIF measurement technique, 
good agreement was found between all sets of calculated 
properties. The spread in the measured values was less for 
a series than for the overall test period. 

Effect of body geometry- The bow shock wave strength 
and properties downstream are defined using the Knudsen 
number (Kn). Kn = X/A is defined as the ratio of the mean 
free path of a gas ( X) to the stand-off distance (A) relative 
to a hemisphere. For frozen flow, such as used in the 
SCFC code, the Knudsen number can be approximated 
from the free-stream Mach number and Reynolds 
number 2. One expression derived by Pai (ref. 28) which 
is applicable from continuum through the transitional 
flow regions is: 

Kn = M oc /(^ 2 ) 1/2 ( 8 ) 


where 

9*2 =(PooU 00 /P2) d H 

The Knudsen number for a flat-faced configuration must 
be corrected to account for the difference in bow shock 
wave stand-off distance between it and a hemisphere. 
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The ratio of bow shock wave stand-off distances for a 
flat-faced configuration and hemisphere plotted as a 
function of the density ratio (p 2 / pi ) across the bow 
shock is illustrated in sketch 3. 
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Sketch 3 


These values were calculated by Katzen and Kaattari, 
assuming equilibrium gas properties (ref. 29). The sketch 
shows that for weak bow shock waves (p2^Pl >0.16) 
the stand-off distance for a flat-faced configuration is 
roughly 2.5 times greater than for a hemisphere with 
the same base diameter. For strong bow shock waves 
(p 2 /pi <0.16) the stand-off distance for the flat-faced 
configuration can increase to a value of greater than 5.0. 

Based on the assumption of frozen chemistry, the 
Rankine-Hugoniot expression was used to calculate the 
density ratio across the bow shock wave in the SCFC 
code (ref. 26). 

p2/P~=(Yf + 0 M i / (Yf-0 M « +2 (9) 

where the free-stream specific heat ratio in the code was 
defined by the frozen state of the gas (ref. 30). 

Yf ~ c p v = ^ y\ c vi 

The density ratio calculated using equation (9) for the 
arc -jet test conditions (table 2) suggested that the bow 
shock wave was relatively weak with a density ratio near 
p 2 /pi = 0.2. Therefore, the Knudsen number for the test 
models was reduced by a factor of 2.6 (sketch 3) from the 
calculated values for a hemisphere to account for the 
difference in their stand-off distances. 

The Knudsen numbers calculated for the flat-faced test 
configurations (cylinder and cone) during the arc -jet tests 
are represented by the square and triangular symbols, 
respectively, in figure 10, which shows that the flow 
environment for the flat-faced cylinder was well within 
the viscous flow region and that the environment for the 


blunt cone was within the merged layer flow region. 
Therefore, Goulard’s theory provides an accurate flow 
simulation for the calculation of the heat transfer rate to 
the test configurations during these arc -jet tests. 

Included in figure 10 is the typical operating envelope 
for both configurations in the AHF using nozzle exit 
diameters from 30.48 to 78.7 cm and enthalpies from 
14.5 to 27.0 MJ/kg. 

Differences in the velocity gradient at the stagnation point 
of the flat-faced test configurations and a hemisphere 
must also be accounted for in Goulard’s theory (eq. (6)). 

Stagnation point heating- In facilities with flows that 
are relatively frozen, the ratio of the velocity gradient at 
the stagnation point of a flat-faced cylinder and hemi- 
sphere has been shown to be directly proportional to the 
ratio of their heat fluxes (refs. 31 and 32): 


(du e / ds) p 
(du e /ds) H 


= (qF /< iH ) 2 


(10) 


where, from Newtonian theory, the velocity gradient for a 
hemisphere is: 

(du e /ds) H =l/R H ^2T(p w -p 00 )/p eo 

Over the operating envelope of the AHF, the relationship 
of the heat flux to a flat-faced cylinder and hemisphere 
with the same base radius was found to be linear (fig. 1 1 ). 
These data show that a linear relationship is valid for 
both cold-wall and hot-wall heat flux measurements. 
Therefore, the velocity gradient for the blunt cone was 
calculated by substituting the slope from this linear 
relationship into equation (6) and accounting for the exact 
location of the sonic point on the cylinder and blunt cone. 
However, the linear relationship in Figure 1 1 would not 
necessarily be valid for all test facilities because of 
possible differences in the state of the gas in the free 
stream and shock layer. 

Finally, calculating the atom recombination coefficients 
for both nitrogen and oxygen required the following basic 
assumptions: (1) a first-order reaction occurs on the 
surface of the material and (2) the accommodation 
coefficient for the material is equal to 1 .0 and NO 
reactions are negligible. 

With the first assumption, the following well known 
expression results: 

Yj =k wi / A /(*T w /2.7iM i ) (II) 


Relative coefficient for air- The SCFC code also 
calculates a relative reaction rate constant (k w ) for each 
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materia] using only air test data and the following semi- 
empirical relationship developed by Rosner (ref. 33): 


2/3 


k w “ (poo U 00 S th )/p w ■ L 


x [(9wo 9 min Vfa max A wo)] 


( 12 ) 


The parameters qmin an( * qmax were obtained from 
Goulard’s theory by setting 0 = 0 and 0 = 1.0, respec- 
tively. The parameter q 0 is the same value used in 
equation (6) (measured radiated heat flux plus heat 
conduction into the model). The recombination 
coefficient (y^) is calculated using equation (11). 

SCFC code validation- In order to substantiate the 
predicted values for the free-stream and shock layer 
properties using the SCFC code, they were compared 
with predicted values using GASP. 

Predicted properties for test cases 2, 3, and 4 are 
compared in figures 12-15. The GASP computation 
supports the frozen flow assumptions used in the SCFC 
code. Predicted free-stream specie concentrations, Mach 
number, and density ratio across the bow shock wave 
calculated from the SCFC code for the nitrogen test cases 
show good agreement with predicted values obtained 
using GASP. 

For both air test cases the free-stream Mach numbers 
calculated using the SCFC code were within 10% of those 
calculated using GASP (figs. 13 and 15). The differences 
in Mach number are attributed to the slight variations in 
the calculated free-stream specific heat ratios (stream 
chemistry) between the two solutions. Therefore, the 
SCFC code coupled with LIF measurements can be used 
to reduce the cost, and the time required for determining 
atom recombination coefficients for candidate SSTO TPS 
is shortened from hours to minutes. 

In addition to calculating the recombination coefficients, 
the code calculates the total hemispherical emittance (e t h) 
using inputs of measured heat flux and surface tempera- 
ture obtained during the arc -jet tests and the well known 
equilibrium reradiation equation: 

£th = 9o ) (1*^) 


Gas properties in the code are obtained from the 
Aerotherm Chemical Equilibrium (ACE) code (ref. 34) 
and using Gupta’s thermodynamic properties (ref. 35). 


Test Results 

Thermochemical Stability 

The effect of high-energy hypersonic flow on the stability 
of the candidate TPS was studied using X-ray fluores- 


cence analyses, mass loss, and emittance measurements 
(figs 16-23). These data indicate that the rigid fibrous 
insulations and carbon composites coated with high 
viscosity glass systems remained relatively stable during 
arc -jet exposure. Carbon composites, with low viscosity 
glass coatings, and preoxidized metallics showed the loss 
of volatile species such as sodium, boron, and/or chrome 
and iron oxides after arc-jet exposure (figs. 16 and 17). 

Removal of these elements is reflected in the mass loss 
data shown in figure 18. For example, the percent mass 
loss for the C-CAT composite sample increased with 
increased surface temperature (figs. 18(a) and 18(b)). The 
apparent weight gain of the LVP sample after an initial 
weight loss appears to be the result of the formation of 
ZrO during the arc -jet exposures (fig. 18(a)). Finally, the 
loss of chrome and iron oxides from the surface of the 
metallic samples had little effect on their mass loss 
histories (fig. 18(c)). The small mass loss experienced by 
these samples is believed to be the result of the initial 
oxides being converted into other forms by the presence 
of the atomic oxygen in the hypersonic air stream. 
Therefore, the metallic samples tested in this study must 
be classified as active TPS. 

Another active TPS is SIRCA, a silicone-impregnated 
reusable ceramic ablator which dissipates the heat at high 
heat flux by mass removal. At low heat flux, where 
ablation does not occur, the heat is dissipated through 
equilibrium reradiation after the formation of a char. The 
char, predominantly a silicon-oxy-carbide (Si*0\ Cy), is 
created through endothermic chemical reactions during 
the arc -jet test. Therefore, its mechanical and surface 
thermal properties (emittance, catalytic efficiency, etc.) 
are not well defined, but depend directly on the specific 
aerothermodynamic heating environment to which they 
are exposed. 

Thermochemical stability of SIRCA- 15F in a 
convectively heated environment can best be illustrated 
by plotting the mass loss rate relative to surface tempera- 
ture (fig. 18). The data used to develop the normalized 
mass loss rates plotted in this figure were obtained from 
reference 17. These earlier tests were also conducted in 
the AHF with similar stagnation point conditions (fully 
dissociated oxygen and partially dissociated nitrogen 
species). Figure 18 shows that SIRCA- 15F starts to ablate 
at a surface temperature near 1 800 K (qHW = 52 W/cm^). 
Included in the figure is a plot of the values for the B^ 
function for SIRCA- 15F. Arrhenius expressions for both 
the mass loss rate and B' w function were obtained from 
curve fits to these data. These expressions, given below, 
can be used in heat shield design calculations. 
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Mass loss rate: 


m/A =(1.0E-7)p°' 8 e a0081ITw (14) 

B^, function: 

B' w = I.6E-8 e 0 00898T w (15) 

Figure 18(b) shows that the function for SIRCA-15F 
is independent of pressure; however, the range of 
pressures was limited during the tests (0.03 to 0. 13 atm). 

Emittance 

Thermochemical instability (loss of volatile species from 
the surface of the TPS) can affect not only the physical 
stability (mass loss, morphology, etc.) but the emittance 
of the TPS (Figs. 20-23). These data show that arc-jet 
exposure of the samples using basically a glass matrix 
affected the spectral emittance throughout the short 
wavelength region (0.3 to 5.0 ji) (fig. 20) and those using 
a basic silicon carbide matrix affected the spectra at the 
longer wavelengths (where the absorption band is 
located) (Fig. 21 ). The absorption band for the coatings 
using a silicon carbide system occurred near 12 microns. 
For the nickel-based systems, the spectral hemispherical 
emittance increased over the entire wavelength range 
because of further oxidization of its surface during arc -jet 
exposure. For the iron-based system (MA-956), the 
emittance was not greatly affected by arc -jet exposure 
(Fig. 23). 

These changes in the spectral reflectance between pre- 
and posttest samples, due to the instability of the surface 
chemistry, have a direct effect on the total hemispherical 
emittance. 

Calculated total hemispherical emittance from posttest 
samples using measured reflectance data is compared 
with values calculated using arc -jet data and the SCFC 
code (figs. 24-26). The values obtained from the 
reflectance data assume that the surface of the sample was 
opaque and they did not include the effect of temperature 
on the spectral hemispherical emittance. In general, the 
two sets of values agreed well for most samples with 
coatings using a basic glass matrix. The values calculated 
from arc -jet data for samples with SiC coating systems or 
oxidized metallic surfaces were lower than those obtained 
from the room temperature spectral reflectance data. 

In fact, for the SiC/SiC and ORCC systems the values 
obtained from the arc -jet data agreed better with the 
pretest calculated values. The spectral reflectance data 
(Figs. 21(a) and 21(b)) show that an increased total 
hemispherical emittance should have resulted after 
arc -jet exposure because of the major effect it had on 
the absorption band for both samples. However, the 


temperature effect on the reflectance data is unknown and 
could have a major effect on the result. One exception 
was the Russian SiC coating system. For this sample, 
both sets of emittance calculations were in good agree- 
ment with each other and, in fact, the trend in the values 
from the arc -jet data support the pre- and posttest results 
obtained from the reflectance data measured at room 
temperature. To fully understand the results of the total 
hemispherical emittance calculations for these TPS 
systems requires further study that includes the effects of 
temperature on the spectral reflectance data. 

The total hemispherical emittance of both S1RCA and 
PM- 1 000 was dramatically changed by arc -jet exposure 
in air (Figs. 24(c) and 26(c)). In both cases the endo- 
thermic reactions on the surface of these samples greatly 
affected the spectral hemispherical emittance, thereby 
resulting in large variations in their total hemispherical 
emittance. SIRCA, which initially has very low emittance 
at temperatures above 700 K, increased to values near 
0.9; and PM- 1000, which initially has very low values at 
temperatures below 1000 K, increased to values near 0.7 
over the entire wavelength range after arc -jet exposure. 
However, for a flight case, the total hemispherical 
emittance of either TPS would not necessarily be the 
same values as obtained in the arc -jet test environment. 
Their emittance could be anywhere between the pre- and 
posttest extremes, depending on the flight trajectory and 
location on the vehicle. 

Total hemispherical emittance for oxidized Inconel 617 
and MA-956 showed little change after arc -jet exposure at 
temperatures above 1000 K. Values calculated using the 
arc -jet data were in close agreement with or slightly lower 
than those calculated from the reflectance data. Expres- 
sions were developed such that at the lower surface 
temperatures they predict values that are close to those 
calculated using room temperature reflectance data, and at 
the higher temperatures they are close to those calculated 
using arc-jet data. These expressions, illustrating the 
temperature dependence of each coated system, are 
summarized in Appendix A. 

Surface Catalytic Efficiency 

Calculated coefFicients for both nitrogen and oxygen 
atom recombination, obtained from data taken in the arc- 
jet and side-arm reactor, are plotted in figures 27-3 1 . 
These values are represented by Arrhenius expressions, 

Yi = ae^^*), which are used in various reacting 
boundary layer or Navier-Stokes solutions to evaluate 
the effect of a material’s catalytic efficiency on surface 
heating (refs. 24 and 36). The slope of the expression (E) 
is representative of the activation energy relative to 
surface reactions. 
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Expressions for RCG, PCC, C-9, TUFI, TABI, and 
preoxidized Inconel 617 were derived from data taken in 
both arc -jet and side-arm-reactor facilities. Expressions 
for the remaining samples were calculated using data 
taken only from arc -jet tests. Coefficients for the candi- 
date TPS can be grouped into three material classes: 
insulators (glass matrix), semiconductors, and metals 
(figs. 27-31). 

For candidate TPS using a basic glass matrix, the 
variations in the coefficients with temperature are similar 
to RCG (fig. 27(a)). Because of viscosity, both nitrogen 
and oxygen coefficients for RCG increased and then 
decreased as the surface temperature approached its upper 
use limit. Similar trends were found in the coefficients for 
grey C-9, C-CAT, and the LVP systems (figs. 28(b), 

29(c), and 29(d)). The characteristic reversal in the 
coefficients for these candidate glass systems occurred at 
temperatures much lower than for RCG. Insufficient data 
were collected on ORCC, precluding complete definition 
of the coefficients for the ORCC system. Coefficients for 
ORCC calculated using arc -jet data suggest that this 
unique characteristic occurs at a surface temperature 
below 1250 K (fig. 29(b)). Further study is required to 
fully define the Arrhenius expressions for ORCC over its 
entire range of use temperatures. PCC (a higher viscosity 
glass system than RCG) coated Nextel 440 was unique 
in that it exhibited the characteristic reversal in the 
coefficients for oxygen but not for nitrogen (fig. 28(c)). 

SiC/SiC, TABI, SIRCA char, Russian SiC, and TUFI are 
grouped as semiconductor materials. Expressions for 
these systems have similar slopes for the oxygen and 
nitrogen coefficients. Note that the apparent activation 
energy (slope) for the oxygen atom recombination 
coefficients is positive and very high (>12,000) for all 
five samples. However, the activation energy for the 
nitrogen atom recombination coefficients for these 
systems is negative and varies over a wide range of 
values. The three systems with porosity (TUFI, SIRCA, 
and TABI) had higher catalytic efficiency than those with 
nonporous surfaces (figs. 27(a), 27(c), and 28(a)). SIRCA 
had the highest value for the oxygen coefficient (y 0 = 0.46 
at 1367 K) and TUFI the highest value for the nitrogen 
coefficient (Yn =0.33 at 1811 K) of these systems. At 
this time it is not clear whether the increased catalytic 
efficiency of all three samples (TUFI, SIRCA, and TABI) 
was only the result of surface porosity or was affected by 
surface chemistry. 

Finally, the calculated values for Inconel 617 are plotted 
in figure 31. Arrhenius expressions used to fit both 
oxygen and nitrogen values have negative slopes and 
increase continuously with increasing temperature. Two 
Arrhenius expressions were required in order to define the 


nitrogen coefficients, but only one to define the oxygen 
coefficients for this preoxidized metallic. Maximum 
values of the coefficient occurred at 1500 K (yn = 0.8 
and Yo = 0.19). 

Arrhenius expressions for all candidate TPS, along with 
the temperature limits for their application, are 
summarized in Appendix B. 

Coefficients derived from Rosner’s expression (eq. (12)) 
are shown in figures 32 and 33. The air coefficients for 
the candidate TPS can also be divided into three groups in 
which their temperature dependences are: 

1 . Not influenced more by the recombination of either 
oxygen or nitrogen atoms. These TPS include 
coatings using basic glass matrices (RCG, C-CAT, 
C-9, and ORCC). 

2. More strongly influenced by the recombination of 
nitrogen than of oxygen atoms (air coefficient 
increases with increased surface temperature). These 
TPS include basically systems that have high silicon 
carbide contained on their surface (TABI, SiC/SiC, 
and SIRCA) and include Inconel 617. 

3. More strongly influenced by the recombination of 
oxygen than of nitrogen atoms (air coefficient 
decreases with increased surface temperature). These 
TPS include TUFI, PCC, and the Russian SiC 
samples. 

Arrhenius expressions for the air coefficients are included 
for each system in Appendix B. Unlike the individual 
coefficients for both nitrogen and oxygen atom recombi- 
nation on an advance TPS, the air coefficients reflect their 
general catalytic efficiency. They are very useful in 
screening candidate TPS materials for application on 
advanced vehicles. However, the individual coefficients 
are required to properly size the TPS for a flight vehicle. 

Surface Heating Distribution 

Arc-Jet Case 

To further validate the SCFC code for use in calculating 
atom recombination coefficients from arc-jet data, surface 
properties obtained from it for RCG (Appendices A 
and B) were used in GASP to calculate the surface 
temperature distributions over a 5 deg blunt cone for two 
arc -jet test cases (table 1, condition 2, for both nitrogen 
and air). The predicted values are compared with 
measured data taken during the tests (fig. 34). The 
predictions agree within 20 K of the measured data 
obtained from the model during the arc-jet tests in both 


11 



nitrogen and air. This is well within the accuracy of the 
temperature measurement. 

Flight Case 

To illustrate the importance of the surface properties of 
individual candidate materials on the overall performance 
of the TPS during Earth entry on a spacecraft, the heating 
distribution along a generic SSTO vehicle’s centerline 
was calculated. The SSTO configuration is simulated 
using a blunt-nosed 33 deg cone with nose radius of 
1 .55 m. The surface of the nose cap is covered with an 
RCG coating, and the fuselage with SiC TABI blankets. 

The heating distribution was calculated using the reacting 
boundary layer code (BLIMPK) which was modified to 
include the wall energy balance for atom recombination 
and the surface emittance effect on the heating (ref. 36). 
These solutions were made for a trajectory point near 
where maximum heating would occur on the SSTO 
winged body at Moo = 22 and at an altitude of 7 1 .5 km 
(fig. 35). 

First, the heating distribution along the centerline of the 
vehicle was calculated assuming equilibrium boundary 
layer flow. The small increase in temperature at the 
junction formed by the two materials is the result of the 
difference in their total hemispherical emittance. Second, 
a BLIMPK solution of the heating distribution along the 
centerline of the vehicle’s midfuselage was calculated 
assuming that both the nose cap and fuselage were 
covered with a borosilicate glass. This solution showed 
that the surface temperatures were much lower than the 
solution assuming equilibrium boundary layer flow. 
Finally, a BLIMPK solution was obtained for the heating 
along the centerline of the vehicle assuming that an RCG 
glass nose cap and silicon carbide TABI covered the 
fuselage. This solution showed that the use of the two 
different materials would result in a large discontinuity 
in the temperature distribution at the junction where the 
glass-coated nose cap and blankets meet. At this junction, 
the surface temperature rapidly increased from a value 
for the RCG coated surface to a value above the one 
predicted assuming equilibrium boundary layer flow. 
Beyond the junction, the temperatures decrease toward 
values predicted assuming an equilibrium boundary layer 
over the surface of the vehicle. This large discontinuity in 
the heating distribution can have a major effect on the 
performance of the vehicle’s TPS as well as its structure 
during an Earth entry and must be accounted for during 
the design. 


Conclusions 

A facility (LASER complex) was described in which 
atom recombination coefficients for proposed TPS can be 
calculated over a wide range of temperatures. Theories for 
calculating the coefficients in both arc -jet and side-arm- 
reactor facilities were described. 

1 . LIF diagnostic techniques were required to accurately 
define the state of the gas during each test in the arc 
jet and specie concentrations in the side-arm reactor. 

2. SCFC code, coupled with LIF measurements, 
provides a very fast, accurate method for determining 
the atom recombination coefficients for candidate 
TPS for future spacecraft. 

3. Accurate predictions of the atom recombination 
coefficients of a candidate TPS require data from 
both side-arm-reactor and arc -jet facilities. 

4. Predicted stagnation point surface temperature 
using atom recombination coefficients for RCG from 
the present SCFC code in GASP agreed well with 
measured values taken from the blunt cone during 
both nitrogen and air tests. 

Surface properties of 16 candidate TPS were defined 
using data and theories developed during this study. 

These TPS include rigid fibrous insulations, reinforced 
carbon carbon, blankets, and metals. 

1 . Test data show that SIRCA and the metals must be 
classified as active TPS. Their thermochemical 
instability is reflected in optical properties and mass 
loss data obtained during the arc -jet tests. 

2. In general, the total hemispherical emittance obtained 
from arc-jet data agreed well with values obtained 
from room temperature-measured reflectance data 
for TPS with a basic glass matrix, but overpredicted 
the values for those with SiC or oxidized surfaces. 

3. Atom recombination coefficients for the candidate 
TPS can be summarized in three groups in which 
their temperature dependences are (1 ) not influenced 
more by either oxygen or nitrogen atom recombi- 
nation (glass systems), (2) more strongly influenced 
by nitrogen atom recombination (SIRCA, SiC/SiC, 
and TABI), and (3) more strongly influenced by 
oxygen atom recombination (TUFI, PCC, and 
Russian SiC). 


12 



Appendix A 

Total Hemispherical Emittance 


Rigid Fibrous Insulations: 

RCG 

Carbon Composite systems: 

SiC/SiC 

e th = 0-89 


e, h = 0.667 + 2.28E -4Tw - 6.8E-8Tw 2 


TUFI 

ORCC 

£ ih = 0.89 


e )h = 0.72 1 +1 .69E -4Tw - 5E -8Tw 2 


**** SIRCA 

C-CAT 

e th = 0.89 (arc -jet exposure) 

e th = 0.846- 1.32E-5Tw 


LVP 


e (h = 0.896 - 1 ,2E -5Tw - 1 ,8E -8Tw 2 


Blanket Systems: 

SiC Coated Russian ACC 

TABI 

e th = 0.766- 1.3 IE -4Tw - 3.2E -8Tw 2 

e, h = 0.836 + 8.27E -5Tw - 1 . 1 1 6E -7Tw 2 

C-9 Coated Nextel 440 


e, h = 0.873 - 1.26E -4Tw + 2.9E - 7Tw 2 -2E - lOTw 3 

Metals: 

PCC Coated Nextel 440 

Oxidized Inconel 617 

e th = 0.875 - 1 .06E-4Tw + 1 ,78E-7Tw 2 - 1 89E - 1 1 Tw 3 

e th = 0.626 - 7.7E -5Tw - 1 .32E -8Tw 2 


MA 956 


e lh = 0. 176 - 3.4E -4Tw - 7.2E -7Tw 2 

PM4000 

e th = 0.757 - 1 .3E -4Tw - 1 ,0E -7Tw 2 
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Appendix B 

Surface Catalytic Coefficients 

Arrhenius expressions for the candidate TPS are listed below. 


Rigid Fibrous Insulations: 


(Tw < 764K) 

RCG 


y N =7E - 6e l443/Tw 

Nitrogen: 


Oxygen: 

(Tw >1575 K) 


(Tw > 1324K) 

y N = 6.2E -6 e*2too/T w 

(RN-1) 

y Q = 1.0E - 6e 14440/Tw 

(905K < Tw <1575 K) 


(685K < Tw < 1324K) 

y N = 10 £-1036 (VT w 

(RN-2) 

y 0 = IOr 6900/Tw 

(465K < Tw < 905K) 


(Tw < 685K) 

y N = 2.0E -5e l500/T w 

(RN-3) 

y 0 = 2.9E - 4e 264/Tw 

(Tw < 465 K) 


Air: 

y N = 5.0E-4 

(RN-4) 

( 1255K < Tw <1644K) 

Oxygen: 


y A = 1 .0E - 4 e 9930/ Tw 

(Tw >1617 K.) 



y 0 = 39E-9e 214l0/Tw 

(RO-1) 

SIRCA 

(978 <Tw <1617 K) 


Nitrogen: 

y G = 5.2 e _8835/Tw 

(RO-2) 

(1255K<Tw< 1644K) 

(502 K < Tw <978 K) 


y N = I.27E -3 E -i6702rrw 

y G = 1.6E -4e 1326/Tw 

(RO-3) 

Oxygen: 

(Tw < 502K) 


(1255K <Tw < 1644K) 

y Q = 5.0E -3e' 400/Tw 

(RO-4) 

y 0 = 7.2E -7e l8283/T w 

Air: 


Air: 

(Tw >1653 K) 


( 1255K < Tw < 1644K) 

y A = 1.5 X 10’ 5 e l2000/T 

(RA-2) 

y A = l.6e' 5250/Tw 

( 1210 K<Tw <1653 K) 



y A = 3.2t" 8285/Tw 

(RA-1) 




Blankets Insulations: 

TUFI 


TABI 

Nitrogen: 


Nitrogen: 

(Tw > 764K) 


(Tw > 1302K) 

y N = 215e- ll730/T ^ 

(TUN-1) 

y N = 580e l4l49/Tw 


(TUN-2) 

(TUO-1) 

(TUO-2) 

(TUO-3) 

(TUA-1) 


(SN- 1) 
(SO-2) 
(SA-3) 


(TAN-1) 
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PCC Coated Nextel 440 


(300K < Tw < 1302) 


Y n = 0.025e* 1006/Tw 

(TAN-2) 

Nitrogen: 


Oxygen: 


(300K < Tw < 1644K) 


(Tw < 1333K) 


Y n = 0.025e' 46(VTw 

(PCN- 1) 

Y 0 = 1.0E -6e I7080/Tw 

(TAO-1 ) 

Oxygen: 


(470K < Tw < 1333K) 


(Tw > 13I6K) 


Y 0 = 5.0e" 3490/Tw 

(TAO-2) 

Y 0 = 6.0E -6e i3040/Tw 

(PCO-1) 

(Tw < 470K) 


(870K < Tw < 1316K) 


Y o = 0.12e- |743/Tw 

(TAO-3) 

Y 0 = I5e' 632VTw 

(PCO-2) 

Air: 


(Tw < 870K) 


(1255K < Tw < 1644K) 


Y 0 = 0.02e‘ 440/Tw 

(PCO-3) 

Y a = 0.45e' 2582/Tw 

(TAA- l) 

Air: 




(1255K<Tw< 1644K) 


Grey C-9 Coated Nextel 440 


Y a = 2.0E -3e 4024/Tw 

(PCA-1) 

Nitrogen: 




(Tw > 132 IK) 




Y n = 1.2E -7e 17034/Tw 

(C9N-1) 

Carbon Composite Systems: 


(930K < Tw < 132 IK) 


SiC/SiC 


Y n = 1.0e" 40I3Tw 

(C9N-2) 

Nitrogen: 


(300K < Tw < 930K) 


(1255K<Tw< 1644K) 


Y n = 0.02e' 340/Tw 

(C9N-3) 

Y n = 10e' 8888/Tw 

(SC N-l) 

Oxygen: 


Oxygen: 


(Tw > 1312K) 


(1255K < Tw < 1644K) 


Y 0 = 2.3E -3e 4877/Tw 

(C90-1 ) 

Y 0 = 4.0E -6e i2348/Tw 

(SCO-1) 

(952K < Tw < 1312K) 


Air: 


Y 0 = 1 .84e' 3900/Tw 

(C90-2) 

(1255K < Tw < 1644K) 


(300K < Tw < 952K) 


Y a = 0.26£- 2895/Tw 

(SCA-1) 

Y 0 = 0.04e" 253/Tw 

(C90-3) 



Air: 


ORCC Coated ACC 


(1255K < Tw < I644K) 


Nitrogen: 


Y a = 3.0E -6e i3540/t w 

(C9A-1 

(1255K<Tw< I644K) 



y N = 1.0E-4e 6623/T w (ORN-1 ) 

Oxygen: 

(1255K<Tw< I644K) 

y 0 = 1 ,0e -7e 15764/Tw (ORO-1) 
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Air: 


Air: 


(1255K <Tw <1644K) 


(Tw > 1567K) 


y A = I.OE -5e 10l65/Tw 

(ORA- 1) 

y A = 3.0E - 3 £ 2653/Tw 

(LPA-1) 



<1255K<TW< 1567K) 


C-CAT 


y A = 0.08e' 25l8/Tw 

(LPA-2) 

Nitrogen: 




(Tw > 1475K) 


SiC Coated Carbon (Russian) 


y N = 6.2E -6 e I2100/Tw 

(CTN-1) 

Nitrogen: 


(l255K<Tw< 1475K) 


(1255K<Tw< 1644K) 


y N = 4.0e' 7625/Tw 

(CTN-2) 

y N = 0.074E -6 e' 236I/Tw 

(RUN-1) 

Oxygen: 


Oxygen: 


(Tw > 1359K) 


<1255K<TW< I644K) 


y 0 = 5.0E -8 e i8023/Tw 

(CTO-1) 

y Q = 4.2E -8E 175337r w 

(RUO-1) 

<1255K<TW< I359K) 


Air: 


y D = 13.5 e _8350/Tw 

(CTO-2) 

<1255K<TW< 1644K) 


Air: 


y A = 8.0E -4e 5040/Tw 

(RUA-l) 

(Tw > 1504K) 




y A = 2.0E -7e l7750yTw 

(CTA-1 ) 



<1255K<TW< 1504K) 


Metal Systems: 


y A = 0.5 e' 44IO/Tw 

(CTA-2) 

Inconel 617 (Preoxidized) 




Nitrogen: 


LVP Coated ACC 


(Tw > 758K) 


Nitrogen: 


y N = 22.0 e‘ 4970/Tw 

(ICN-I) 

(Tw > 1529K) 


(300K < Tw < 758K) 


y N = 1.5E -Se 10080 ^* 

(LPN-l) 

y N = 0.1 1 e -953/Tw 

(ICN- 2) 

<1255K<TW< 1529K) 


Oxygen: 


y N = 0.06 e' 2605/Tw 

(LPN-2) 

(300K < Tw < 1644K) 


Oxygen: 


y 0 = 0.55 e ,6IWTw 

(ICO-1) 

(Tw > 1499K) 


Air: 


y 0 = 2.5E -7e 17533/T w 

(LPO-1) 

<1255K<TW< 1478K) 


(1255K<Tw< 1499K) 


y A = 1.7E +4E-l 52757r w 

(ICA-1) 

y 0 = 7.5e’ 8283/Tw 

(LPO-2) 
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Table 1 . Candidate TPS 


Coating 

Substrate 

Sample 

Test models 

Supplier 

l.RCG 

FRCI-12 

Blunt cone/ 
split tube 

Blunt cone/ 
split tube 

NASA/ARC 

2. TUFI 

AETB-12 

Blunt cone/ 
split tube 

Blunt cone/ 
split tube 

NASA/ARC 

3. SIRCA 

FRCI-15 

Blunt cone 

Blunt cone 

NASA/ARC 

4. None 

TABI 

Disk/cloth 

F-F cylinder/ 
tube liner 

NASA/ARC 

5. Grey C-9 

Nextel 440 

Disk/cloth 

F-F cylinder/ 
tube liner 

Rockwell 
International Corp. 

6.PCC 

Nextel 440 

Disk/cloth 

F-F cylinder/ 
tube liner 

NASA/ARC 

7. None 

Nextel 440 

Cloth 

Tube liner 

NASA/ARC 

8. C-CAT 

Reinforced 

carbon-carbon 

Disk 

F-F cylinder 

NASA/LRC 

9. LVP 

Reinforced 

carbon-carbon 

Disk 

F-F cylinder 

NASA/LRC 

10. ORCC (RT-4) 

Reinforced 
carbon -carbon 

Disk/split tube 

F-F cylinder 

ROHR 

1 1 . None 

SiC/SiC 

Disk 

F-F cylinder 

McDonnell Douglas 
Aerospace 

12. Russian 

Reinforced 

carbon-carbon 

Disk 

F-F cylinder 

McDonnell Douglas 
Aerospace 

13. None 

Inconel 617 

Disk/sheet 

F-F cylinder/ 
tube liner 

McDonnell Douglas 
Aerospace 

14. Oxide 

Inconel 617 

Disk/sheet 

F-F cylinder/ 
tube liner 

McDonnell Douglas 
Aerospace 

15. Oxide 

PM- 1000 
(nickel-based) 

Disk/sheet 

F-F cylinder 

ROHR 

16. Oxide 

MA-956 

(iron-based) 

Disk/sheet 

F-F cylinder 

ROHR 
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Table 2. Typical arc -jet operating conditions 


Test case 

Power, MW 

Mass flow, gm/s 

Pressure, atm 

Enthalpy, MJ/kg 





Bulk/Effective 

Nitrogen 

I 

1.12 

23.6 

0.68 

12.4/13.2 

2 

2.8 

45.4 

1.36 

18.0/20.3 

3 

4.2 

69.4 

2.04 

19.6/20.9 

4 

5.44 

101.2 

2.72 

19.1/19.8 

Air 

1 

1.04 

23.6 

0.68 

12.1/13.9 

2 

2.52 

45.4 

1.36 

17.5/18.8 

3 

5.0 

101.2 

2.72 

17.1/19.8 

4 

7.2 

157.9 

4.08 

16.9/18.6 


20 



21 


Sketch of complex 






(a) Side-arm reactor 


SWING-IN 

CALORIMETER 



(b) Aerothermodynamic heating facility 


Figure 2. Schematics of facility equipment 






(a) Conical sample 



(b) Cylindrical sample holder 


Figure 3. Test models. 
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PMT Signal 




(a) Low catalytic efficiency sample, long tube 




(b) Medium to high catalytic efficiency sample, composite tube section 
Figure 4. Typical fluorescence (PMT) signal distribution along diffusion tube in side-arm reactor. 
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Nozzle Solution 


Shock Layer/Model Solution 


A*,M« 


(a) Typical axisymmetric grid for Navier-Stokes computation using GASP 




(b) Frozen Mach number flow computation using effective area ratio 
Figure 5. Control volume for free-stream and shock layer simulations . 
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5600 



Frozen Mach Number, Mf 

(a) Free-stream velocity in nitrogen (H^ =19 MJ/kg, P ch = 4.08 atm) 



Frozen Mach Number, Mf 


(b) Atomic nitrogen in free stream of arc-jet flow 
Figure 6. Predicted free-stream properties relative to frozen Mach number. 
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c 

CJ 
c J 

e 

o 

U 

E 

o 


> 


QJ 

a: 



Energy, AH, Joules/cc 


faj Relative nitrogen atom concentration 



Chamber Pressure, atm 


(b) Free-stream temperature 


Figure 9. Variation of stream properties in the AHF during surface catalysis studies 


Knudsen Number, X/A 


10 F 


I Merged Layer 


!4 


0.1 


0.01 


0.001 


Viscous Layer 


I- Vorticitx Interaction 


IT 


C Continuum 



Squares - Blunt Cone 
Triangles- Flat-Faced Cylinder 


P ch = 0.68 atm 
Heo = 27.9 MJ/kg 


P ch =4.06 atm 
Heo = 14 MJ/kg 


Operating envelope for AHF 


10 


100 1000 10000 

Effective Area Ratio, A^A* 


Figure 10. Flow field in front of arc-jet test models during exposures to air. 



Heat flux, W/cm 2 
(Flat-faced cylinder) 


Figure 1 1. Correlation of stagnation point heat flux between hemisphere and flat-faced cylinder. 

























Arbitrary units 



(a) Posttest photograph MA-956 



Post test surface chemistry (center) 

(b) Posttest photograph PM-1000 

Figure 17. Pre - and posttest photographs and surface chemistry of metallic TPS. 
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Percent mass loss, AM/Mj xlOO 









00 

© 



(a) Mass loss rate 



Wall Temperature, K 
(b) Dimensionless mass transferrate 


Figure 19. Mass loss characteristics of SI RCA during arc-jet exposure in air. 
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Spectral Hemispherical Emittance 


Pre-test 


Post test 



Q I .. .1 u . I i . . I. tL li ii I ii iLU il *-«J ll *l lt.Il-1 

0 2 4 6 8 10 12 14 16 !8 20 0 2 4 6 8 10 12 I 4 16 18 20 


(a) C-9 coated Nextel 440 



0 2 4 6 8 10 12 I 4 16 18 20 0 2 4 6 8 10 12 I 4 16 18 20 


(b) PCC coated Nextel 440 



Wavelength, microns 

(C) SIRCA 

Figure 20. Effect of arc-jet exposure on spectral emittance of fibrous TPS. 
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Spectral Hemispherical Emittance, 


Pre-test 











Total Hemispherical Emittance, 





Total Hemispherical Emittance, e, 


Dashed lines - calculated from room temperature spectral reflectance data 



(a) SiC/SiC (b) ORCC 



(c) C-CAT (d) LVP 



0 500 1000 1500 2000 


Surface Temperature, K 

(e) Russian ACC 


Figure 25 . Total hemispherical emittance of carbon composite TPS. 
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Total Hemispherical Emittance, e, 


Dashed lines - Calculated from room temperature spectral reflectance data 



0 500 1000 1500 2000 



0 500 1000 1500 2000 


(b) MA-956 



0 500 1000 1500 2000 


Surface Temperature, K 
(c) PM -1000 

Figure 26. Total hemispherical emittance of metallic TPS 
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Atom Recombination Coefficiently 


Nitrogen 



0.0004 0.0012 0.002 0.0028 0.0036 


Oxygen 



0.0004 0.0012 0.002 0.0028 0.0036 


(a) RCG 



0.0004 0.0012 0.002 0.0028 0.0036 0.0004 0.0012 0.002 0.0028 0.0036 


(b) TUFi 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 


1/Temperature, K' 1 
(c) SIRCA 

Figure 27. Atom recombination coefficients for rigid fibrous TPS. 
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Atom Recombination Coefficient, 


Nitrogen 

Circles- arc-jet data 



Oxygen 

Triangles - Side-arm-reactor data 



( a ) TAB I 


1 B* 


: i Yn = 1.2E-7E 
:\ / y N = l.oe- 4013 ^ 0 

\ Ym =oj 


0.1 

0.01 

0.001 


0.0001 


17034/T(o 


y N =0.02£- 34O/To) 


0.1 

0.01 

0.001 


i .y 0 = 2.3E“3E 4877/Tc) 

| V * / Yo = I.84e- 3900/Tco 

. * * 


Li-LLUJUlJ 


lilt 1.L-I.1-LJ 


0.0001 


\ 


1 y 0 = 0.04e' 


• * ■ ■ 1 ' ■ 1 ■ ««■»*«»* i-J 


-253/T(o 


1 1 i 


0.0004 0.0012 0.002 0.0028 0.0036 


0.0004 0.0012 0.002 0.0028 0.0036 


(b) Grey C-9 coated Nextel 440 




1/Temperature, K* 1 
(c) PCC coated Nextel 440 

Figure 28. Atom recombination coefficients for blanket TPS. 
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Atom Recombination Coefficient, 


Nitrogen 


Oxygen 


1 n- 


0.1 t 


0.01 t 


o.ooi L 


1 r 


X 


y N = ioe^ 888/Tt ° 


0.1 

0.01 

0.001 


/ 


7o = 4.0E-6e 


12348/Tto 


0.0001 i ■■■ * » — i * i i i— i — * * i 0 0001 

0.0004 0.0008 0.0012 0.0016 0.0004 


1 - 1 - 1 


0.0008 


0.0012 


0.0016 


(a) SiC/SiC composite 



(b) ORCC 



1/Temperature, K* 1 
(c) C-CAT coated ACC 

Figure 29. Atom recombination coefficients for advanced carbon composites. 
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Atom Recombination Coefficient, 


Nitrogen 


Oxygen 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 

(a) LVP coated ACC 



1/Temperature, K' 1 


(b) Russian ACC 


Figure 30. Atom recombination coefficients for advanced carbon composites (concluded). 
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Atom Recombination Coefficient 



0.0004 0.0012 0.002 0.0028 0.0036 

1/Temperature, K' 1 
(a) Nitrogen 



1/Temperature, K- 1 

(b) Oxygen 

Figure 31. Atom recombination coefficients for Inconel 617. 
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Atom Recomination Coefficient 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 


(a) RCG (b) TUFI 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 


(c) I SIRCA (d) TABI 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 


(e) C-9 (Nextel 440) (>) PCC (Nextel 440) 

1/Temperature, K-* 

Figure 32. Air recombination coefficients for both rigid and blanket TPS. 
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Atom Recombination Coefficients 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 

(a) SiC/SiC (b) ORCC 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 



0.0004 0.0008 0.0012 0.0016 0.0004 0.0008 0.0012 0.0016 


(e) SiC (Russian) (f) Inconel 61 7 

1/Temperature, K' 1 

Figure 33. Air recombination coefficients for both carbon composites and metals . 
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Surface temperature (K) 


Measured surface temperature 



(a) Nitrogen flow 


<$> Measured surface temperature 



(b) Air flow 


Figure 34. Predicted and measured surface temperature on blunt cone during test case 2. 
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